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Abstract. Combinatorial optimization is a fertile testing ground for statistical physics methods developed 
in the context of disordered systems, allowing one to confront theoretical mean field predictions with actual 
properties of finite dimensional systems. Our focus here is on minimum matching problems, because they 
are computationally tractable while both frustrated and disordered. We first study a mean field model 
taking the link lengths between points to be independent random variables. For this model we find perfect 
agreement with the results of a replica calculation, and give a conjecture. Then we study the case where 
the points to be matched are placed at random in a d-dimensional Euclidean space. Using the mean field 
model as an approximation to the Euclidean case, we show numerically that the mean field predictions are 
very accurate even at low dimension, and that the error due to the approximation is 0(l/d 2 ). Furthermore, 
it is possible to improve upon this approximation by including the effects of Euclidean correlations among 
k link lengths. Using k = 3 (3- link correlations such as the triangle inequality), the resulting errors in the 
energy density are already less than 0.5% at d > 2. However, we argue that the dimensional dependence of 
the Euclidean model's energy density is non-perturbative, i.e., it is beyond all orders in k of the expansion 
in fc-link correlations. 

PACS. 75.10.Nr Spin-glass and other random models - 02.60.Pn Numerical optimization 



1 Introduction 

1.1 Background 

The study of disordered and frustrated systems, and 
in particular spin-glasses, has long been a major issue 
in condensed matter physics (for reviews see [0,|^,§]). 
Most efforts have been based on replicas, in part because 
that method has led to the exact solution of the 

Shcrrington-Kirkpatrick (SK) model ||. However, since 
the SK model is of infinite range, it is not clear @,|].f| 
how relevant its solution is for understanding finite dimen- 
sional spin-glasses such as the Edwards- Anderson (EA) 
model |lCj] . The application of the replica formalism to fi- 
nite dimensional systems, on the other hand, is hampered 
by two major difficulties jll],|l^]. First, in the saddle point 
equations, finite connectivities lead to an infinite number 
of order parameters: one has to deal with order parameters 
Qap-y... having an arbitrarily large number of indices ]l3| . 
In contrast, the SK model requires only the order parame- 
ter q a p, with two indices. The second difficulty comes from 
the Euclidean nature of space: The metric structure intro- 
duces constraints on the possible values of the quenched 
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disorder variables between points % and j. In the infinite 
range model, these variables are independent, but in the 
short range models, the values allowed depend on the dis- 
tance between points i and j, thereby introducing strong 
correlations. These correlations make the replica analysis 
much more difficult. 



In this article we study matching problems; these 
are disordered and frustrated short range models aris- 
ing in combinatorial optimization. They are simpler than 
spin-glasses and the difficulties just stressed have to a 
large extent been overcome. In particular, for the mini- 
mum matching problem based on independent random link 
lengths, Mezard and Parisi (M & P) have worked out the 
order parameters <7 aj g 7 ... exactly JhJ] and have introduced 
a link correlation expansion pa] to take into account cor- 
relations among link variables. Using extensive numerical 
simulations of the ground states, we study the following: 
(i) the probability distribution of ground state energies in 
the disorder ensembles; (ii) the distribution of "local" en- 
ergies in the ground state configurations; (iii) the validity 
of their replica approach in the case without correlations; 
(iv) the accuracy of this random-link length model as a 
mean-field approximation to the Euclidean model; (v) the 
accuracy of the link correlation expansion. 
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1.2 The models 

Consider A points (N even), and a set of links with lengths 
lij = Iji connecting points i and j, = 1,... , A). We 
call an instance a specification of these link lengths. A 
matching of the points is a set of bonds in which each of 
the A points is the extremity of one and only one bond. In 
effect this is a dimcrization: points are linked in pairs. The 
length or cost of a matching is then the sum of the lengths 
of its bonds. (Since the points will often belong to a Eu- 
clidean space, we will use the term length rather than cost 
hereafter.) The minimum matching problem (MMP) jl6| 
can be defined as the problem of finding a matching of 
minimum length, given the hj's. A variant of this problem 
is the minimum bipartite matching problem (MBMP) : 
we now have two sets of A points instead of one (A is no 
longer necessarily even) and a bipartite matching consists 
of a dimerization where each pair contains one point from 
each set. The length and the minimization of the bipartite 
matching are defined as for the MMP. From an algorithmic 
point of view, both problems belong to the class P, mean- 
ing that the exact solution of any instance of size A can 
be found in a time at most polynomial in A. In fact, the 
MMP can be solved in a time of 0(A 3 ) and the MBMP 
in a time of 0(N 2 In A). This property is important for a 
numerical study since it allows extensive exact computa- 
tions. We have implemented the 0(NE\n(E)) matching 
algorithm as exposed in [[l7[ where E is the number of 
edges in the graph. The resulting program solves a typical 
A = 200 random point instance in less than a second on 
a Dec-Alpha machine. The total computation time spent 
generating the data summarized in the tables of this paper 
amounts to about 100 days of machine time. 

These combinatorial optimization problems may be 
mapped onto physical systems, where each matching is 
one state of the system with an energy equal to its length. 
Then the minimum matching problem is equivalent to 
finding the ground state of that system. The physical sys- 
tems built in this way are frustrated since in general not 
all points may be matched with their nearest neighbor. 
The thermodynamics of these systems can be studied by 
introducing Boltzmann weights for the states as was pro- 
posed by Kirkpatrick et al. |Q. Since this article focuses 
on ground states, however, we restrict our discussion to 
the zero temperature properties of these systems. 

In any disordered system, be it a spin-glass or a match- 
ing problem, one is not particularly interested in the prop- 
erties of a given instance of the problem. More relevant 
physically are typical properties, or averages over an en- 
semble of instances. The Z^'s then become quenched ran- 
dom variables and one speaks of the stochastic M(B)MP, 
problems that are both disordered and frustrated. More- 
over, we are interested in the infinite "volume" limit, 
meaning in this case the limit A — > oo. These systems can 
then be studied using the replica or the cavity method (for 
a review see 0) as developed by M & P 

Let us now describe the different ensembles of lij's 
we consider. The first is the Euclidean ensemble (for the 
MMP as well as for the MBMP) where we have A (or 
2A) random and independent points chosen uniformly 



in a d-dimensional volume (e.g., a unit hypercube), and 
the lengths Zy are given by the Euclidean distances be- 
tween the points. This is a short-range model in that the 
points lie in a Euclidean space and only the first few near- 
est neighbors are relevant for the minimum matching. (In 
the formulation based on spins Q , the coupling between 
spins decreases exponentially with distance.) Using a spin- 
glass analogy, this Euclidean model is like the EA model 
(both models are short range, leading to important Eu- 
clidean correlations in the disorder variables) . Not surpris- 
ingly, the Euclidean MMP — like the EA model — is not 
solvable analytically. This suggests that one should con- 
sider a different ensemble for the disorder variables in or- 
der to render the problem more tractable [^0| . Indeed, this 
may be accomplished by taking the lengths iy to be in- 
dependent, identically distributed random variables. The 
corresponding model is called the "random-link" model. 
(Note that the points do not lie in a metric space and 
the ijj's no longer have correlations; in particular the tri- 
angle inequality does not hold). Pursuing our spin-glass 
analogy, the random-link MMP is like the model of Viana 
and Bray |H| for spin-glasses. Both models are "infinite 
dimensional" in the sense that there is no underlying ge- 
ometry and thus there are no Euclidean correlations. Fur- 
thermore, in both cases, the effective connectivity at each 
site stays fixed as the size of the system grows. For the 
Viana and Bray model, this is enforced by having the num- 
ber of non-zero couplings to a spin be size independent; 
for the random-link MMP, this occurs because only the 
first few nearest neighbors of a site effectively contribute 
to the minimum matching. These models thus interpolate 
between the infinite range / infinite connectivity case and 
the finite range / finite connectivity case. Thus, since it is 
expected that the Viana and Bray model provides a bet- 
ter approximation than the SK model to the EA model, 
the fact that its analogue here (the random-link MMP) is 
exactly solvable is of major interest. 

A connection between the Euclidean and the random- 
link models was first given by M & P ]l9| : they pointed 
out that the one- and two-link distributions could be made 
identical in both problems, leading to a particular fam- 
ily of random-link models parametrized by d, where d 
is the spatial dimension for the corresponding Euclidean 
model. With such a choice, both models have the same 
Cayley tree approximation. One may then consider the 
random-link MMP to be a mean-field model for the Eu- 
clidean MMP in which correlations between link lengths 
(the quenched disorder variables) have been neglected. 
What we call from now on the "random-link approxima- 
tion" 
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consists of using the thermodynamic functions 
of the random-link models as estimators for those of the 
Euclidean models. 



1.3 Outline 

This paper expands upon previous work [ p3| , and provides 
an in depth study of ground state properties (disorder in- 
duced distributions of the energy, length of dimers, dimen- 
sional dependence, etc...) in the MMP and the MBMP. 
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The outline is as follows. In section |], we examine the 
random- link models. First, we discuss self- averaging and 
large TV scaling properties of the ground state energy. Sec- 
ond, we recall results derived assuming replica symmetry 
and check them numerically. Third, we propose a conjec- 
ture for the frequency with which a point connects to its 
fcth nearest neighbor in the ground state of the M(B)MP. 
In section |L we study the Euclidean models. First, we dis- 
cuss self-averaging of the ground state energy. Second, we 
quantify numerically the precision of the random-link ap- 
proximation. Third, we consider the Euclidean corrections 
to this approximation in the case of the MMP. Finally, 
section ^ discusses our results and provides an outlook on 
possible generalizations to other systems. 

In Appendix |A|, we present an TV-independent upper 
bound for the length of the random- link MMP when d = 1 . 
In Appendix |b[ we prove a self-averaging property for the 
Euclidean MBMP in dimensions greater or equal to 3. 



2 Random-link models 

In order to use the random-link models as an approxima- 
tion for the Euclidean ones, we must, as previously men- 
tioned, set the random-link distribution to match the one- 
link Euclidean distribution. In the Euclidean model, one 
possible approach would be to take the large TV limit at 
fixed density of points, in which case the volume would 
scale linearly with TV. However, for historical reasons, the 
standard practice is to take TV — > oo in a fixed volume. 
These two pictures are equivalent, and are mapped onto 
one another by rescaling all the lengths by a factor N 1 ^. 
As a consequence, in the units we shall use in the rest of 
this paper, the mean length between neighboring points 
scales as 7V~ 1 / d . 

Consider two points i and j randomly chosen in the 
d-dimensional unit hypercube. In the absence of edge ef- 
fects, the distribution of their distance lij is given by 
Pd(kj = r) = dB d r d -\ where B d = n d / 2 /{d/2)\ is 
the volume of the d-dimensional unit ball. We thus take 
from now on pd{l) as the individual distribution for the 
link lengths in the random-link model. The random-link 
M(B)MP models are then described by a single parameter 
d. Two comments are in order. First, p d is not normalized, 
and so must be cut off; this can be done arbitrarily as the 
large TV scaling of the minimum matching depends only on 
the behavior of pd(l) at small I. Second, neglecting edge ef- 
fects, any two lengths are uncorrelated in Euclidean space, 
so the prescription just given also matches the two-link 
Euclidean distribution. 



scale as N~ x / d , so L^f M is expected to scale as N 1 ~ 1 / d . 
For instance, consider the length di between point i and 
its nearest neighbor. It is easy to show that in the large 
TV limit, (di) ~ 73i(d)/A 1 / d where 



D 1 (d) = (l/d)\B d 



l/d 



(1) 



is the rescaled average nearest neighbor length at large 
TV. (In the units where the density of points is equal to 
one, Di(d) is exactly the mean nearest neighbor length.) 
In the case of the minimum matching length, it can be 
proven |24| that L^ M jN x ^ 1 l d becomes peaked around 
its mean value as TV — > oo. A stronger property, called 
self-averaging, would be that L^/TV 1-1 /'' tends with 
probability one to a (non-random, TV-independent) con- 
stant P M L ^(d), as TV — ► oo. Although this second prop- 
erty has not yet been proven, it is strongly supported 
by previous numerical studi es 2j| as well as the simula- 
tions described in paragraph 2.3. Moreover, bounds can be 



found for (3 M L M (d) (see Appendix |a|) that reinforce this hy- 
pothesis. In the following, we assume the existence of the 
i^mmM' The analogous discussion applies to MBMP as 



well, so we will assume (3 



RL 

MBM 



[d) 



Urn L 



RL 

MBM 



/N 



l-l/d 



It is of interest to understand the distribution of L^ M 
in the large TV limit. In particular, one may wonder 
whether the self-averaging property comes from some kind 
of central limit theorem which would lead to a Gaussian 
(normal) limit distribution for L^ M . The central limit 
theorem (CLT) states that for a sum S of TV indepen- 
dent identically distributed random variables the relative 
standard deviation a = ^JVar(S) / (S) and the skewness 
s = ((S - (S)) 3 )/<t 3 decrease as 1/y/N. Although in fact 
the terms entering the sum L^fj^ are correlated, if the 
correlations are not too strong, we can expect the same 
CLT-type scaling to hold for Lf^ M . We have computed a 
and s numerically for Lf;f M and find a behavior that is in 
excellent agreement with the expected CLT scaling laws 
(see Table El). (Since the finite size effects and the statis- 
tical noise are significant, we have not extrapolated our 
data to the TV — > oo limit.) 

Similarly, we have computed a and s for the MBMP, 
and find nearly identical results (see again Table 0). 
Clearly, the large TV behavior follows the CLT scaling laws 
for both the MMP and the MBMP. Furthermore, our data 
suggest that the amplitudes associated with these scalings 
are equal in the two models. We conjecture that the two 
are equal for any of the connected mo ment s, and provide 
theoretical support for this in section |]^. We will also 
derive a theoretical estimate for the values of a and s in 



section 2.3 



2.1 The large TV limit 

We denote by L^ M the length (or energy) of the mini- 
mum matching in the MMP; L^ M is a random variable 
depending on the instance (i.e., on the i^'s). It is a sum 
of TV/2 terms, each of which is typically the length be- 
tween near neighbors. We have seen that these lengths 



2.2 Survey of analytical results 

The random-link MMP was first solved by M & P [juj 
using the replica method, with a replica symmetric ansatz. 
They verified p^j the stability of the replica symmetric 
solution (at least for d = 1), suggesting that most likely 
the ansatz is exact. They also confirmed their results by 



4 



J. Houdayer et al.: Comparing Mean Field and Euclidean Matching Problems 



Table 1. Numerical measurements of the relative standard 
deviation a and skewness s of the distribution of L^m an d 
LkFbm at d = 1 and d = 2. 





MMP d 


= 1 


MMP d 


= 2 


N 


osJW]2 s 


y/N/2 


s 


s/W[2 


50 


0.784 


1.39 


0.421 


0.33 


100 


0.798 


1.50 


0.423 


0.32 


200 


0.798 


1.39 


0.427 


0.36 


400 


0.808 


1.56 


0.424 


0.24 




MBMP d 


= 1 


MBMP d 


= 2 


N 




sVN 




sVW 


25 


0.780 


1.29 


0.419 


0.25 


50 


0.795 


1.36 


0.423 


0.29 


100 


0.802 


1.36 


0.427 


0.35 


200 


0.811 


1.57 


0.428 


0.45 



the cavity method p9| . In our notation their results may 
be written 



dDi(d) 



2(l/d)\ J_ 

where Gd(x) satisfies the integral equation 

r + oc 



G d (x)e- G ^ x) dx, (2) 



G d (x) = d (x + y) 



' d{v) d% 



(3) 



and Di(d) is as defined in equation [l| Equation [}] can 
be solved at d = 1, leading to G\{x) = ln(l + e x ) and 
(3 M L M {1) = tt 2 /24. Furthermore, M & P calculated the 
distribution Pd(l) of the rescaled bond length N x / d lij in 
the minimum matching in the limit N — > oo, and found 



P d (l)=dl 



d-1 



dG 



dx 



d {x)e~ G ^- G ^ l - x Ux. (4) 



In the case d = 1, they found Pi (I) — 2(21 — 
e~ 21 sinh(22))/ sinh 2 (2Z). Finally, since replica symmetry 
is not broken, one expects the mean length to have a 1/7V 
expansion: 



(L 



RL \ 
MMI 



oRL 

Pmm 



(d) 1 



A(d) B(d) 



N 



N 2 



(5) 



M & P have calculated E6| the first subleading term at 
d = 1 and have found « —0.13. 

Similar calculations were performed for the MBMP 
(Orland ||7), M & P Q[P}||2§]). In our units, they find: 

P^BM(d) = W R AAd)- (6) 

At d = 1, they obtain (3 M L BM (1) = tt 2 /12 and A(l) m 
-1.61. 

The MBMP and the MMP are very closely related. 
This is best seen within the convention we have used 
where the MMP matches N points, the MBMP matches 



2N points and both models have the same individual link 
length distribution One then sees that the saddle point 
equations for the partition functions (in the limit of large 
N) become identical in the two models. Thus Pd(l) is the 
same in the two models, and in fact any reasonable ob- 
servable will be the same in both models at large N. This 
remarkable property has apparently gone unnoticed so far. 
(Given this property and our conventions, the factor 2 in 
equation |^ simply follows from the fact that a bipartite 
matching has twice as many bonds as a simple matching.) 
A consequence of this correspondence is that the moments 
of L M L M and L M L BM should be the same at large N. We 
are thus able to provide theoretical support for the con- 
jecture, given in the previous section, and based on our 
numerical results. 

It is also of interest to study the limit of large d. We 
have derived a l/d expansion of (3 M L M (d). One way to do 
this is to set G d (x) = Gd(x = x/d + 1/2) and then write 
Gd(x) as a power series in l/d. From this we find |2~4| 



Pmm (d) 



Di(d) 



1_7 vr 2 /12 + 7 2 /2- 7 /l 

d d 2 U 3 



where 7 = 0.5772 ... is Euler's constant. 



(7) 



2.3 Numerical verifications and a new conjecture 

Brunetti et al. p5| have used numerical simulations of the 
random-link models to confirm the predictions of (3 RL to 
the level of 0.2% for the MMP and 0.7% for the MBMP 
at d = 1 and d = 2. They have also checked the 0(1 /N) 
corrections to (L RL ) at d = 1 and find relatively good 
agreement with the theory. 

To obtain further confirmation we have estimated 
f3 M L M {d) and Pmbm id) numerically for 1 < d < 10, and 
have found accordance with the replica symmetric pre- 
dictions to the level of 0.03% for the MMP and to the 
level of 0.1% for the MBMP. In order to reduce the sta- 
tistical fluctuations and get quantitative errors on our es- 
timates for j3, we use the following procedure. First we 
compute the ensemble average (L R f M )/N 1 ~ 1 ' d using a 
variance reduction trick |2^,^8|. (The values of N used 
here and for the other calculations of (3 given later on are 
N = 50, 70, 100, 150, 200, 260, and 400; we performed av- 
erages over 5 ■ 10 5 instances at N — 50 down to 7, 500 
instances for N = 400.) Then, to get the large N limit, we 
fit our data using a 1/N series (truncated after the second 
order) as indicated by the theory (equation ||). The fits 
are good, with y 2 values confirming the finite-size scaling 
law. The statistical error bar on (i M L M (d) is then obtained 
by the standard method |g{J, whereby fixing the fit's lead- 
ing coefficient at (3 ± a makes x 2 increase by one from its 
minimum value. Our results are summarized in Table |. 

We have also checked the prediction for the distribu- 
tion Pd(l) of the bond lengths in the minimum (bipartite) 
matching. We find such good agreement with theory that 
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Table 2. Comparison of theoretical and numerical values of 
Pfj L M (d) and /SmbmM- Numbers in parentheses represent the 
statistical error bar on the last digit(s). 







MMP 


MBMP 


d 






j3num 1 2 


1 


0.411234 


0.41142(12) 


0.41134(9) 


2 


0.322580 


0.32257(5) 


0.32262(4) 


3 


0.326839 


0.32684(4) 


0.32691(4) 


4 


0.343227 


0.34323(3) 


0.34327(4) 


5 


0.362175 


0.36210(3) 


0.36222(3) 


6 


0.381417 


0.38143(3) 


0.38143(3) 


7 


0.400277 


0.40026(5) 


0.4002(5) 


8 


0.418548 


0.41852(5) 


0.4180(5) 


9 


0.436185 


0.43612(5) 


0.4362(5) 


10 


0.453200 


0.45310(5) 


0.4531(4) 



the numerical data is indistinguishable from the replica 
predictions (on a figure one would not be able to tell the 
two curves apart). 

We will now use Pd(l) to obtain an estimate of the 
quantities computed numerically in Table [I]. Neglecting 
the correlations among the bond lengths in the optimal 
matching, we can use equation |] to predict a and s as 
a function of A. At d = 1, we find VyjN/2 m 0.73 and 
sy/N/2 m 1.85 for the MMP. For the MBMP, the same 
formulae apply if A/2 is replaced by N. These theoretical 
predictions (to be compared with Table [j]) are about 10% 
too small for a and about 30% too large for s, showing 
that correlations cannot be neglected, but are nevertheless 
relatively small. 

Another interesting quantity is the mean fraction of 
points connected to their fcth nearest neighbor in the op- 
timal matching. Call this fraction Pd{k). In view of our 
numerical data at d = 1 (see Table ||), we conjecture for 
both the MMP and the MBMP that in the limit N -> oo, 



Pi(fc) =2 



-Us 



(8) 



In addition the A-dependence of these fractions seems to 
be linear in 1/A as one could expect from equation [| It 
seems likely that our conjecture may be confirmed by us- 
ing the replica method. There may also be analogous re- 
lations in higher dimensions, but unfortunately we have 
not found any convincing formulae. At best, our data 
are approximately fitted by stretched exponentials (see 
also ]28|). This same kind of behavior also arises in other 
combinatorial optimization problems such as the traveling 
salesman problem |50|. It is also interesting to note that 
another conjecture has been proposed for a random-link 
MBMP by Parisi His conjectured relation gives the 
mean length of the matching at all values of A, not just 
in the limit A — > oo. 



3.1 The large A limit 

Let L MM be the length of the minimum matching in 
the Euclidean MMP. Hereafter, we take the points to 
be distributed randomly in the d-dimensional unit hyper- 
cube. Following the same argument as for the random- 
link MMP, one expects Lf IM to scale as N 1 ~ 1 / d . In fact, 
it has been proven J3^] that the Euclidean MMP has the 
self-averaging property in any dimension, so Lf IM /N 1 ~ 1 ' (i 
tends to a constant (3 MM W as A — > oo with probability 
one. 

For the Euclidean MBMP, the situation is more com- 
plex. There are two sets of points, so local density differ- 
ences have a large effect. In particular, the 7V 1-1 / d scal- 
ing law is not valid for d < 2. At d = 1, it is easy to 
see that the optimum corresponds to matching the points 
left to right. Then a quick estimate shows that L MBM 
scales as VA instead of as A . Furthermore L MBM / VA 
does not become peaked, so there is no self-averaging. At 
d = 2, the situation is more interesting: L MBM scales [[53) 
as %/Aln A (instead of as VA for Lff BM ). The question 
of self-averaging has not yet been settled but numerical 
simulations indicate that the property does hold, and we 
will assume this is the case hereafter. At higher dimen- 
sions, we have proven self-averaging (see Appendix [Bh, so 



lim 



N- 



^MBM 



/ N l-l/d ex j sts for 



the quantity (3f [BM (d) 
d>3. 

As with the random- link models, one may wonder 
whether the central limit theo rem is at work. Following 
what was presented in section 2.1, it is natural to inves- 



tigate the limiting distribution of the optimum length. 
It is convenient in the numerical study to avoid bound- 
ary effects; to do so, we work in the unit hypercube with 
periodic boundary conditions. Using numerical measure- 
ments of the moments of the distributions we find that 
the MMP obeys the CLT scaling laws (see Table |). This 
is not surprising: for the Euclidean MMP, we can divide 
the hypercube into subvolumes. Then the length of the 
minimum matching is close to the sum of the minimum 
matching length in these subvolumes, and the CLT scal- 
ings should hold. On the other hand, the CLT scalings do 
not hold for the MBMP, and furthermore the limit distri- 
bution is not Gaussian. In particular, we find that s does 
not tend to zero, rather it grows with A. The CLT argu- 
ment used for the MMP does not apply to this problem 
because the subvolumes just mentioned will not in general 
contain equal number of points from each set. Neverthe- 
less, as the dimension increases, the density fluctuations 
decrease, explaining why at fixed A the MBMP values of 
a and s get closer to the MMP ones as d increases. (For a 
formal application of this argument, see Appendix [Bl) 



3.2 Numerical results and the random-link 
approximation 



3 Euclidean models 



To date, little has been done to compute the ground-state 
energy densities (3f IM (d) and 0f IBM (d) . The best esti- 
mates prior to our recent work Q were |54|,^8| (3f 4M (2) sa 
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Table 3. Numerical results of 10 5 (pd(k) — 2 k ) for the random-link MMP and MBMP in the case d = 1. Numbers in parentheses 
represent the statistical error on the last digit (s). 



MMP 

N = 50 N = 100 N = 200 



MBMP 

N = 50 N = 100 JV = 200 



1 


422(22) 


181(16) 


86(16) 


1185(16) 


608(11) 


305(11) 


2 


174(19) 


96(14) 


47(14) 


113(14) 


47(10) 


17(10) 


3 


-26(15) 


-12(10) 


1(10) 


-238(10) 


-108(7) 


-52(7) 


4 


-81(11) 


-30(8) 


-18(8) 


-289(7) 


-157(5) 


-77(5) 


5 


-112(8) 


-49(5) 


-31(5) 


-248(5) 


-123(4) 


-56(4) 



Table 4. Numerical results for the relative standard deviation a and skewness s of the distribution of Lf /M and Lff BM at 
d = 2, 3 and 4. 





MMP d 


= 2 


MMP d 


= 3 


MMP d 


= 4 


N 


s 




s 




osJW[2 s 




50 


0.302 


-0.90 


0.244 


-0.65 


0.206 


-0.68 


100 


0.299 


-0.93 


0.244 


-0.66 


0.205 


-0.64 


200 


0.290 


-0.87 


0.243 


-0.52 


0.204 


-0.52 


400 


0.295 


-0.85 


0.243 


-0.88 


0.203 


-0.45 




MBMP d 


= 2 


MBMP d 


= 3 


MBMP d 


= 4 


N 




sVN 


aVW 


sVN 




sVW 


25 


0.603 


2.61 


0.350 


1.16 


0.258 


0.42 


50 


0.744 


4.20 


0.387 


1.94 


0.274 


1.12 


100 


0.938 


6.49 


0.431 


3.26 


0.289 


1.88 


200 


1.197 


9.97 


0.481 


4.71 


0.305 


2.53 



Table 6. Comparison of MBMP ground state energies for the 
two models: Euclidean and random-link (3 < d < 10). The 
rightmost column compares the Euclidean MBMP and MMP 
models. For f3 E (d) the numbers in parentheses are statistical 
errors on the last digit(s). In the last column Ae = (Pmbm ~ 
2/?mm) / 1 Pmbm ■ 



d 


P E {d) 


P RL {d) 


j3 RL_ j3 E 


d 2 Ae 


3 


0.7080(2) 


0.653679 


-7.68% 


0.664 


4 


0.7081(2) 


0.686455 


-3.06% 


0.597 


5 


0.7349(1) 


0.724350 


-1.44% 


0.514 


6 


0.7688(3) 


0.762834 


-0.77% 


0.482 


7 


0.8039(2) 


0.800554 


-0.41% 


0.461 


8 


0.8391(2) 


0.837097 


-0.24% 


0.445 


9 


0.8736(2) 


0.872370 


-0.14% 


0.429 


10 


0.9076(2) 


0.906400 


-0.13% 


0.436 



0.312 and /?mm( 3 ) ~ °- 318 ; for PmbmW, no valid es- 
timates have yet been published. Expanding upon the 
work in |23f| , we now provide very accurate measurements 
of these quantities, using the same procedure as in the 
random-link case. Again we find that the \ 2 values justify 
the use of a truncated 1/N series. 

As previously remarked, our random-link distributions 
were established in order to match the one- and two-link 
distributions of the Euclidean model. So, if the effects 



coming from the Euclidean correlations among three or 
more link lengths are small, then the properties of the 
random-link and Euclidean M(B)MP should be quantita- 
tively close. In fact, replacing /3^ M (ci) and /3ff BM (d) by 
Pfhi W an d Pmbm (d) leads to a very precise approxima- 
tion. As shown in Table ||, P^ L M (d) differs from (3f IM (d) by 
17.8% at d = 1 and by 3.9% at d = 2, and this difference 
decreases quickly as the dimension increases (the qua ntit y 
(3 EC given in the table will be discussed in section |3.3[). 



Likewise for the MBMP, shown in Table |6|: /3$ L BM (d) dif- 



fers from (3i 



MBM\ 

J MBAi(d) by only 7.7% at d = 3. Note that com- 
paring Pmbm($) an< ^ ^mbm(^) at (i < 2 is meaningless 
since the scaling law s ar e different as we mentioned in the 
beginning of section 3.1. Nevertheless, we have also com- 



puted /3f, IBM (2); empirically, we found that using a 1/7V 
fit could not do, but that good \ 2 values were obtained us- 
ing a linear fit in l/\n(N). We find PmbaA 2 ) = 0.340(1). 

Since the random-link approximation was also found to 
be very good for the traveling salesman problem ^] , our 
results suggest that this approximation should be widely 
applicable to link-based optimization problems. Further- 
more, we can understand how the size of the error inherent 
to this approximation decreases as d — » oo. Consider for 
instance the bond occupation probabilities for links con- 
nected to a given site. One expects that these probabilities 
have a large d limit, and that Euclidean correlations intro- 
duce 1/d corrections compared to the random- link case. 
(To be precise, one expects the corrections to be given 
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Table 5. Comparison of MMP ground state energies for the three models: Euclidean, random-link, and random-link including 
3-link Euclidean corrections (1 < d < 10). For f3 E (d) the numbers in parentheses are statistical errors on the last digit. 



d 


nE ! t\ 

0(d) 


p id) 




,0 RL -0 E 
ft H H 


aEC ( j\ 

P (d) 




fi E 




fi E 


1 


0.5 


0.411234 


-17.75% 


-0.178 






2 


0.3104(2) 


0.322580 


3.92% 


0.078 


0.30915 


-0.40% 


3 


0.3172(2) 


0.326839 


3.04% 


0.091 


0.31826 


0.33% 


4 


0.3365(3) 


0.343227 


2.01% 


0.080 


0.33756 


0.30% 


5 


0.3572(2) 


0.362175 


1.39% 


0.070 


0.35818 


0.27% 


6 


0.3777(1) 


0.381417 


0.98% 


0.059 


0.37849 


0.21% 


7 


0.3972(1) 


0.400277 


0.77% 


0.054 


0.39807 


0.22% 


8 


0.4162(1) 


0.418548 


0.56% 


0.045 


0.41685 


0.17% 


9 


0.4341(1) 


0.436185 


0.48% 


0.042 


0.43485 


0.17% 


10 


0.4515(1) 


0.453200 


0.38% 


0.038 


0.45214 


0.14% 



by a 1/d expansion, with the leading term naturally be- 
ing of order 1/d.) In support of this argument, we have 
checked that the relative difference between the Euclidean 
value and the random- link value of Pd(k) is indeed of or- 
der 1/d. We can then see how quantities such as (3(d) 
depend on the link lengths vary when one goes from the 
random-link model to the Euclidean model. At large d, a 
simple calculation shows [^2| that the mean link lengths 
between fcth nearest neighbors have for different values of 
k a relative difference of 0(1 /d) (regardless of whether 
we use the random-link or the Euclidean ensemble) . Since 
the occupation probabilities Pd{k) have Euclidean correc- 
tions of order 1/d, the random- link approximation leads 
to an error of order 1/d 2 for (3(d) and Pd(l). We have 
tested this behavior numerically by considering the quan- 
tity d((3 M L M — MM )//3 MM , an d the data scale as 1/d as 
expected. Thus the random-link approximation gives both 
the leading and subleading 1/d dependence of Pf IM (d). 
Then, from equation ^, we find 



1.2 



Pmm (d) 



DM 



(9) 



Furthermore we have directly confirmed this dependence 
by performing a fit of our (3^j M (d) values and we find 
0.424 ± 0.008 for the coefficient of the 1/d term; this is to 
be compared to the theoretical value 1 — 7 = 0.42278 .... 

Performing the same analysis for the MBMP, we find 
that (3 MBM (d) also satisfies Equation ^ (omitting the fac- 



tor 1/2) or equivalently that (f3 MBM 
0(1/ d 2 ) (see Table |). A direct fit to thcT/d term of 
Equation for Pf, IBM gives 0.42 ± 0.05 in agreement with 
the theoretical value 1 — 7. 

Finally, we have also applied the random-link approxi- 
mation to compare Pd(l) in the random- link and Euclidean 
models and found very good agreement (see Figure |l| for 
the MBMP at d — 3, which has the largest discrepancy). 
The agreement becomes better as the dimension increases, 
with an error of order 1/d 2 . 



0.8 



0.6 



0.4 



0.2 




0.5 1 1.5 



2.5 



3.5 



Fig. 1. Comparison of the random- link (continuous line) 
and Euclidean (points) integrated distributions of the rescaled 
bond lengths in the minimum bipartite matching for the 
MBMP at d = 3 and N = 60. 



3.3 Euclidean corrections to the random-link 
approximation 

The mean field model (i.e., the random- link model) pro- 
vides a very good approximation to the finite dimensional 
Euclidean model. It is of major interest to push the ap- 
proximation further and to derive, for instance, a large 
dimensional expansion. For the MMP, M & P have calcu- 
lated |]l5| a correction to the random-link approximation 
by considering the effects of the three-link Euclidean cor- 
relations. Such correlations arise only when three links 
form a triangle. (For the MBMP, one would have to go 
to four-link correlation effects, and this has not yet been 
attempted). M & P's result is given in terms of a function 
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Gd, but where Gd now satisfies a much more complicated 
integral equation (equation 34 in their paper). From this, 
one obtains a new estimate for (]f tM (d), which we denote 
here by (3ff M (d) (EC stands for Euclidean corrections). 

We solved numerically this modified integral equation 
for Gd and computed (iffL (d) for 2 < d < 10 (see Ta- 
ble ||). Comparing with we see that these esti- 
mates are considerably more accurate than when using 
f3fr[ L M (d). At d — 2, the random-link approximation leads 
to an error of 3.9%; this error is reduced by nearly a factor 
of 10 by including the corrections due to 3-link correla- 
tions. At d = 3, the error is reduced from 3.0% to 0.4%. 
At larger d, the error continues to decrease, though the 
effect is less significant. 

To understand the dependence of these corrections 
on dimension, it is useful to consider how the difference 
PkFm ~ 0mm sca l es with d. This difference is associated 
with a 3-link correction term which gives the probabil- 
ity of finding nearly equilateral triangles as d — > oo. It is 
not difficult to see that this probability goes to zero ex- 
ponentially with d. Thus the 3-link correlations give tiny 
corrections at large d (as confirmed by the numerics), and 
the power series expansion in l/d of (i M c M is identical to 
that of @ M L M . In fact the l/d series of (3ff M is not modi- 
fied if one includes 4, 5, or any finite number of multi-link 
correlations. This follows from the fact that all non-zero 
multi-link correlations come from sets of links forming at 
least one loop, and that fixed sized (iV-independent) loops 
connecting near neighbors become exponentially rare at 
large d. (Of course this behavior depends on having ran- 
domly placed points; the situation is very different on lat- 
tices, for instance, where small loops are important.) The 
main consequence of the nature of these multi-link corre- 
lations is that the M & P fc-link correlation expansion will 
not converge towards (3f 4M (d), and in particular it does 
not even allow one to compute the 0(1/ d 2 ) term in the 
l/d expansion of /3 MM (d). Formally, the Euclidean dimen- 
sional dependence is beyond all orders of the fc-link corre- 
lation expansion. This non-perturbative behavior is quite 
remarkable, and indicates that the large N limit and the 
fc-link correlation expansion do not commute. This prop- 
erty may have its analogue in other disordered systems. 

4 Summary and Discussion 

In this article we have studied two versions of the stochas- 
tic minimum matching and minimum bipartite matching 
problems: the random-point Euclidean ensemble, and its 
correlation-free approximation, the random-link ensem- 
ble. For both ensembles, we have given evidence that the 
ground state energy is self-averaging and obeys (except for 
the Euclidean MBMP) the central limit theorem scaling 
laws. For both ensembles we have performed extensive nu- 
merical simulations in order to measure the ground state 
energy density (3(d) and the distribution Pd(l) of bond 
lengths in the ground-state. For the random-link model 
we have checked to high precision the replica symmetric 
prediction of Mezard and Parisi and find excellent agree- 
ment. Furthermore, we have proposed a new conjecture 



at d = 1 that suggests further analytical calculations. For 
the Euclidean model, we have studied the accuracy of the 
random-link approximation and find the error to be small 
even at low dimensions. For example at d = 2 the error 
introduced by this approximation is 3.9% for (if IM . We 
have also been able to go beyond the random-link approx- 
imation by applying the formalism of Mezard and Parisi 
in order to include 3-link correlations associated with the 
triangle inequality. We find that this improved estimate re- 
duces the error by nearly a factor of 10 at low dimensions. 
In particular, the resulting d = 2 prediction for (3 MM has 
an error of 0.4%. 

The limit of high dimensions is also of major interest. 
Based on our simulations up to d = 10, we have given 
strong evidence that the first two terms of the l/d ex- 
pansion of the random-link model are the same as for the 
Euclidean model but that these two models differ at order 
l/d 2 . Furthermore we have argued that for any fixed k, 
fc-link correlations in Euclidean space only give rise to ex- 
ponentially small contributions in d and thus do not mod- 
ify the l/d expansion. As a consequence, the Euclidean 
dimensional dependence is beyond all orders of the fc-link 
correlation expansion. 

Although our study was performed in the context of 
the MMP, our reasoning applies also to other link-based 
problems associated with random points, and leads us to 
the following picture. For these random-point systems, 
whenever the thermodynamic functions depend only on 
the local properties of the (short) link graph, we expect 
the error in the random-link approximation to be expo- 
nentially small in d. This will always be the case in the 
high temperature phase where correlations are weak. How- 
ever, in the low temperature phase, the correlations may 
be such that the N — * oo limit and the fc-link correla- 
tion expansion do not commute. We expect this to be the 
case in many combinatorial optimization problems such as 
the assignment problem and the traveling salesman prob- 
lem, where fc-link correlations with fc growing in N remain 
important as N — > oo. Arbitrarily large loops matter in 
these systems, and contribute to the thermodynamics at 
order l/d 2 . This change in behavior can be illustrated us- 
ing more physical language by considering a polymer on 
a random-point lattice. Then the arguments given pre- 
viously indicate that the random-link approximation is 
exponentially good in the "dilute" phase, while it leads 
to an error in powers of l/d in the collapsed phase. This 
non-perturbative behavior of the l/d expansion is remi- 
niscent of what occurs in lattice systems such as the Eden 
model, where the l/d expansion does not commute with 
the N — > oo limit |35[ . We would not be surprised if similar 
phenomena occurred in other disordered systems, whether 
they be based on links or spins; the calculation of their l/d 
series would then be particularly difficult. Further insight 
into these issues may be obtained by looking at excited 
states of the M(B)MP to see whether they consist of ar- 
bitrary large loops (i.e., whose sizes diverge as N — > oo). 
Such a study could also allow a direct investigation of pos- 
sible replica symmetry breaking in the Euclidean case. 
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It is tempting to speculate on how all of this might 
carry over to spin-glasses. There, one analogue of the 
random-point MMP is the Edwards- Anderson model with 
nearest neighbor interactions on a <i-dimensional hyper- 
cubic lattice of connectivity 2d. The mean field model 
for spin-glasses is usually taken to be the Sherrington- 
Kirkpatrick infinite connectivity model. However, a more 
appropriate mean field model in this context is that of 
Viana and Bray J2l| in which the connectivity is 2d and 
for which the couplings Jy have the same individual dis- 
tribution as in the Edwards-Anderson model. Just as in 
the matching problem, one thereby obtains a mean field 
model at any given dimension, and the mean field approx- 
imation then consists of using this model to estimate the 
thermodynamic quantities of the d-dimensional Edwards- 
Anderson model. We expect this approach to lead to errors 
of a few percent at low dimensions, and to 0(1/ d) errors 
at large d. Unlike the random-link MMP, the Viana-Bray 
model has not yet been solved analytically, so that the 
mean field values would have to be calculated approxi- 
mately. Nevertheless, we view these mean field models as 
providing a very promising approach to computing quan- 
tities in the Edwards- Anderson model. We hope that the 
potential reward will encourage new attempts to solve the 
Viana and Bray model, and that the challenge of deter- 
mining the first 0(1/ d) correction to "mean field" will be 
taken up. 

We are grateful to D. Dean, W. Krauth, M. Mezard, H. Orland, 
N. Sourlas and J. Vannimenus for comments, and to O. Bohi- 
gas for his constant interest. Special thanks go to A. Percus 
for constructive criticisms. J.H. and J.H.B.dM. acknowledge a 
fellowship from the MENESR, and O.C.M. acknowledges sup- 
port from the Institut Universitaire de France. The Division de 
Physique Theorique is an Unite de Recherche des Universites 
Paris XI et Paris VI associee au CNRS. 



A Bounds for the MMP 

Here we review exact lower and upper bounds for /?mm(^)' 
and provide to our knowledge the first finite upper bound 
at d = 1. First there is a trivial lower bound: each point 
is at best linked to its nearest neighbor, so 

ffikAd)> D 1 (d)/2. (10) 

To get an upper bound, we use another optimization prob- 
lem: the traveling salesman problem (TSP), which consists 
of finding a minimum length tour visiting all the N points. 
Indeed we can obtain a matching by removing every sec- 
ond bond in a tour (note that this argument fails for the 
MBMP), so we get 



&{d) < (3^ P (d)/2, 



(11) 



where PtspW is tne analog of l3^j L M (d) for the TSP. A 
bound for P^g P {d) is already known for d > 2 (see in |20[]). 
This bound comes from the greedy algorithm and is 



tsp (d) 



< 



DM 
1 - 1/d' 



(12) 



In the case d = 1, the greedy construction leads to a bound 
which grows logarithmically in N and thus is useless. We 
have obtained a finite bound for the case d = 1 using a 
different approach. The idea is that a bound is known [|36| 
p7| for the asymmetric TSP (where hj and Iji are indepen- 
dent). Let us call Ay the link lengths of the asymmetric 
TSP. We denote by pAsym(l) the distribution density of 
the Ay, and associate a symmetric TSP to any asymmet- 
ric TSP by setting Zy = min(Ay, Xji). Then, in the limit 
of short link lengths, the distribution in this symmetric 
TSP is p S ym(l) = 2p Asym (l). This gives 

4s7(d) < 2 1 ' d ^-{d). (13) 
And we thus obtain our bound at d = 1. 

B Self-averaging for the Euclidean MBMP 

We shall denote by Lmbm(A"i, . . . , Xjy, Yx, ■ ■ ■ , Yjy) the 
length of the minimum bipartite matching between the 
Xi's and the V^'s. We prove here that for any d > 3, Lmbm 
satisfies a theorem analogous to the one Beardwood, Hal- 
ton and Hammersley have proven for the TSP |58) . Specif- 
ically, we prove the following: 

Let X\, . . . , V/v, ■ ■ ■ and Y\, . . . , Yjsr, ... be two se- 
quences of random points independently and uniformly 
distributed in [0, l] d , where d > 3, and let Ln — 
Lmbm{Xi, . . . , Xpj; Yi, . . . , Y/v). There exists a constant 
PMBM(d) > such that with probability one, 



urn — - — — rr 

N^oc V 1 " 1 / d 



P. 



MBM 



(d). 



(14) 



To begin with, we remark that to prove this theorem, 
it is sufficient to establish that Ln /N 1 ^ 1 / 4 converges in 
mean value to a constant (3m bm (d) . This is a consequence 
of the following lemma ]39|: For any t > 0, 



P 



L 



N 



L 



N 



>t < 2 exp(- 



N l-2/d t 2 



(15) 



The theorem then follows easily from the convergence of 
{Lpf) /N 1-1 / 4 as N — > 00, by applying the Borel-Cantelli 
lemma. We now wish to establish that for d > 3 the 
quantity (Ljv) /V 1_1 / d indeed converges to a constant 
PMBM{d) > 0. To do this, we exploit the subadditivity 
properties of Lmbm (see 40 ). 



First we need to generalize Lmbm to matchings be- 
tween two sets of different cardinalities. We shall define 
Lmbm{X\,... ,X Nl ;Yi,... ,Vv 2 ) by requiring that the 
matchings contain as few unmatched points as possible, 
that is we leave |Vi — V2I points unmatched. 

Suppose the points Xi, . . . X^ 1 , Y\, . . . Ypj 2 belong to 
an arbitrary cube Q whose edges have length a, and di- 
vide Q into disjoint cubes Q p , p = 1, . . . , 2 d by splitting 
each edge in two halves. Construct in each Q p an optimal 
matching in the sense just defined, between the n\ tP points 
Xi and the ri2 lP points Yi in Q p , and denote its length by 
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L p . There will be |roi )P — n 2 , P \ points left unpaired in each 
Q p , so if Lq denotes the length of an optimal matching for 
these points, one has 

2 d 

Lmbm iX\ , . . . , Xp? 1 ; Y\ , . . . , Yjv 2 ) < L p + Lq 

p=i 

2 d ^ 2 d 

|. (16) 

p— 1 " p— i 

Now we apply this to Q = [0, l] d . Let Q Pl ,pi = l,...2 d 
be the cubes obtained in the last subdivision, let Q PlP2 be 
the cubes obtained by splitting in two halves the edges of 
each Q Pl , and so on. By repeating this operation K times, 
we get a subdivision with cubes Q Pl ... PK whose edges are 
of length 1/2 K . Let n l pi PK and n2 :P1 ... PK be respectively 
the number of points Xi and Yi in Q Pl ... PK . Apply ( jl6| ) first 
to the Q Pl ,...p K _ 1 's, then to the Q Pl ... PK _ 2 , s, etc., keeping 
at each step only those points that are still unpaired. It 
is easy to convince oneself that the number of unpaired 
points in each Q Plt ... lPK _ k just after step k is given by 
\ n i,pu— >PK-h ~ n 2,pi,... ,p K _ fc |- After step k = K one ob- 
tains a matching between X\,... , X^ 1 and Y\, . . . , Yiv 2 
where all but \N\ — N 2 \ of the points are matched. One is 
thus led to the following inequality: 

Lmbm(Xi, . . . Xn 1 ;Yi, . . . Yjv 2 ) < ^ L P1 ,,, PK 

P\..-PK 

+ 5Z"2F \n>i, Px ... Pk - n 2 , Pl ...p k \. (17) 

k=l Pi-.-Pk 

We now proceed to derive a subadditivity property for 
the average value of Lmbm- To do this, it is useful to 
consider the case where N± and N 2 are not fixed inte- 
gers but are independent Poisson random variables with 
the same fixed parameter N. For a given k, the num- 
bers ni tPlj ... Ph and n 2 , Pl ,... Pk are then also independent 
Poisson random variables, with parameter N/2 kd . Let 
M(N) = (L MBM (Xi, . . . ,X Nl ;Yi, . . . , Y Na )). It is easy 
to see that 



(L pl ... PK ) = 2- K M(N/2 Kd ) 



(18) 



Moreover, well known properties of Poisson variables allow 
us to write 

/ n\ 1/2 

(\ni, Pl ... Pk - n 2 , Pl ... Pk \) < V2 i . (19) 



By taking mean values in (JIT]), we are thus led to 

K 

M (N) < 2 K( - d - 1 hd(N/2 Kd ) + V2dNj2 2k{d/2 ~ 1] ■ ( 20 ) 



fc=i 



Extending this construction, one may easily prove 
that for 2 K < m < 2 K+1 we have 



A" 



M(N) < m d - 1 M(N/m d ) + 2 d V2dN^2 k{d ' 2 - 1 \ (21) 



Dividing this last inequality by TV 1 1 / d and then replacing 
N by m d N, 



M{m d N) M{N) 



< 



K 



(m d N) 1 -V d ~ N 1 - 1 ^ 7Vi/2-i/d 



2 d V2d 9 -fe(d/2-i) 



(22) 



Standard arguments may now be used to show that 
the ratio M(N)/N 1 ~ 1 / d necessarily converges to a limit 
0MBM(d) as N -► oo. Indeed, let f {t) = M(* d )/t d_1 . One 
verifies at once that f(t) satisfies 



f(mt) < f(t) + C d /t d / 



2-1 



(23) 



for all t > and all integer m; f(t) is continuous, since 
M(N) is a continuous function of TV. So the expression 
fit)^rC d /t d / 2 - 1 is bounded in [1,2] and since [1, oo) is the 
union of the intervals m[l, 2], m > 1, it follows from ( ]23| ) 
that f{t) remains bounded as t — * oo, thus lim* f(t) < oo. 
Now define (3 — lim* /(£). For any e > 0, choose to ^> 1 
and 77 > such that /(t) + C d /t d /' 2 - 1 < f3 + e for t in the 
interval I = [to — r], to + r/]. Since the intervals m/, m > 1 
span a whole interval [A, 00) for A sufficiently large, it 
follows again from ( p3| ) that lim* f(t) < + e. Since e is 
arbitrary, we have lim fit) = /3, hence f(t) — > /? as i — > 
00, from which it follows that limTv^oo M(A r )/A^ 1_1//<i = 
/3. □ 

We have thus shown that for d > 3, 



\Lmbm\X\ 



,%)> 



/3wBw(d)^ 1_1/<i . A^-oo (24) 



when TVi and A2 are independent Poisson variables with 
parameter N . The same result for the mean value (Ljv), 
where N is a fixed integer, then follows easily. Indeed, one 
has the obvious bound 

\Lmbm(Xi,Xn;Yi,Yn) — Lmbm(Xi, Xjy^Yi, Yjv 2 )| 
< Vd(|JVi - AT| + |AT 2 - N\), (25) 

and taking mean values and dividing by N 1 ~ 1 / d , we de- 
duce that 



lim ^rrn = 0mbm{<1). 



(26) 



fe=0 



For further discussion of self-averaging proofs see |32|,|4l| 
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